Coherent vortex structures in quantum turbulence 
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This report addresses an important question discussed by the quantum turbulence community 
during the last decade: do quantized vortices form, in zero-temperature superfluids, coherent struc- 
tures similar to vortex tubes in ordinary, viscous turbulence? So far the evidence provided by 
numerical simulations is that bundles of quantized vortices appear in hnite-temperature superfluids, 
but from the interaction with existing coherent structures in the turbulent (viscous) normal fluid, 
rather than due to the intrinsic superfuid dynamics. In this report we show that, in very intense 
quantum turbulence (whose simulation was made possible by a tree algorithm), the vortex tangle 
contains small coherent vortical structures (bundles of quantized vortices) which arise from the 
Biot-Savart dynamics alone, and which are similar to the coherent structures observed in classical 
viscous turbulence. 



It has been known since the 1980's that homogeneous 
isotropic turbulence contains intermittent worm-like re- 
gions of concentrated vorticity [H-Q- Their role in the 
dynamics of turbulence is not clear Q: there is no con- 
sensus as to whether they are responsible for the main 
properties of turbulence (e.g the celebrated Kolmogorov 
energy spectrum) or only affect the tails of statistical 
distributions and the exponents of high-order structure 
functions. Turbulence is also studied at temperatures 
near absolute zero in superfluid helium [ol-fl^. Here the 
viscosity is zero, and the vorticity is concentrated in thin, 
discrete vortex lines of fixed core radius and circulation 
(whereas vortices in ordinary fluids are viscous and can 
have any size and strength). This 'quantum turbulence' 
shares importa nt prop erties with the ordinary (classical) 
turbulence Il3l. [l4 |. from the same drag force on a 
moving sphere [l5| to the Kolmogorov energy spectrum 
in homogeneous isotropic turbulence |16l421| . Here we 
show that, alongside the Kolmogorov spectrum, quan- 
tized vortices tend to form coherent structures, similar 
to the vortex tubes of ordinary turbulence. Despite its 
relative simplicity, quantum turbulence seems to share 
important features of generic turbulent flows. 

At temperatures T below 0.7 K, thermal excitations 
become ballistic (in that sense that their propagation is 
not affected by any interaction with the background) so 
that liquid ^He is effectively a pure superfluid. Because 
of the quantum-mechanical constraints on the rotation, 
its flow is everywhere potential with the exception of vor- 
tex line singularities of atomic thickness (the core radius 
is approximately ao ~ 0.1 nm), each line carrying one 
quantum of circulation k = h/m — 10~^m^/s, where 
h is Planck's constant and m is the mass of the single 
helium atom. (In liquid '^He-B, which becomes super- 
fluid at much lower temperatures, T < 1 mK, the core 
radius is about 80 nm, and the quantum of circulation is 
K = h/(2m) = 0.662 x lO^^m^/s, where m is the atomic 
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mass of the ^He isotope, whose nuclei are fermions.) 
Quantum turbulence, readily created by stirring the liq- 
uid helium, takes the form of an apparently disordered 
tangle of such vortex lines. The question which we ask 
is whether discrete quantized vortex filaments organize 
themselves and form coherent structures, as in ordinary 
turbulence. 

To model quantum turbulence in the zero-temperature 
limit, we observe that oq is many orders of magnitude 
smaller than the typical distance, £ « 10~*-10~^m, be- 
tween vortex lines in experiments, so it is appropriate to 
describe vortex lines as space curves s(^,t) of infinitesi- 
mal thickness which move according to the Biot-Savart 
law 

where s is the position vector of a point on the vor- 
tex line, t is time, ^ is the arc length measured along 
the vortex line, and the integral extends over the entire 
vortex configuration C Equation ([T|) states classical in- 
compressible Euler dynamics in integral form. However, 
unlike classical Euler vortices, quantized vortices recon- 
nect if they come sufficiently close to each other, as seen 
directly in recent experiments p3| and demonstrated us- 
ing the Gross-Pitaevskii equation for the Bose-Einstein 
condensate To apply equation ^ to superfluid he- 
lium we must therefore add an algorithmic rcconnection 
procedure [25j . 

We perform numerical simulations of reconnecting vor- 
tex lines at T = evolving according equation ([1]) in a 
period cubic box of the size D = 0.75 x 10"'^ m. The 
technique to discretize the vortex lines, the regulariza- 
tion of the Biot-Savart integral, and the vortex recon- 
nection procedure are standard and well described in the 
literature [l^ . Our numerical algorithm, which con- 
trols the discretization, enforces a separation between 
points on the vortex line to remain between 2 x 10~^ m 
and 4 x 10~® m during the evolution; this is therefore 
the numerical resolution of the computation. The time 
integration uses a third-order Adams-Bashforth scheme; 
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the typical time step is At = 4 x 10~^ s. Our vortex 
reconnection procedure guarantees that a small amount 
of vortex length (as a proxy for kinetic energy) is lost at 
each reconnection, in agreement with numerical calcula- 
tions of vortex reconnections performed using the Gross- 
Pitaevskii equation [2BI • The reconnection procedure and 
the enforcement of the minimum distance provide the 
numerical dissipation mechanism which plays the role of 
phonon emission in actual ''He [l^l (or of the Caroli- 
Matricon f28| mechanism in '^He-B). To cope with the 
large number N ^ 10^ of discretization points along the 
vortex lines, we use a tree algorithm [2l|, |2^, [s^l which re- 
duces the CPU time needed to evaluate the Biot-Savart 
integrals from N'^ to A^logA^. The initial condition con- 
sists of randomly oriented straight vortex lines, each car- 
rying one quantum of circulation (however, we checked 
that results do not depend on the initial configuration of 
vortices). We find that the initial vortices interact and 
reconnect, and, after a transient of the order of 0.03 s, 
the vortex line density (vortex length per unit volume) 
saturates at the average value L = 7.9 x 10* m~^. The 
average curvature saturates too. The vortex tangle is 
shown in Figure [TJ By construction, the circulation of 
each vortex line is conserved, and each vortex moves ex- 
clusively due to the collective influences of all the other 
vortex lines; this represents an accurate model of quan- 
tum turbulence. 




FIG. 1: The instantaneous form of the saturated vortex tangle 
at t = 0.06 s. 

One of the most important quantity in theory of tur- 
bulence is the energy spectrum, E{k) defined by 

E^T7 [ ^W\'dV= r E{k)dk, (2) 
^ Jv Jo 

where E is the kinetic energy (per unit mass), V is 
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FIG. 2: The energy spectrum E{k) for the flow of Figure [T] 
The red dashed line shows the classical Kolmogorov scaling 



volume, k — |k|, with k the wave vector. Figure [5] 
shows a spectrum typical of quantum turbulence. The 
Biot-Savart interactions of vortex lines over length scales 
larger than the average inter-vortex distance £ « L^^^^ 
has induced the same Kolmogorov energy spectrum 
E{k) ~ (for k <^ kt = 27r/£) which is observed 

in ordinary turbulence (46j . The Kolmogorov spectrum 
was also observed in experiments with turbulent super- 
fluid helium [TgI . Il7j , and in calculations performed with 
both the vortex filament model [l^ [2l| and the Gross- 
Pitaevskii equation [l^, [l^ . At larger wavenumbers the 
spectrum becomes shallower due to the Kelvin- wave cas- 
cade along individual vortex lines [13, [Sll, [s^l , in agree- 
ment with Refs. [IB, iSSj. The Kelvin cascade (which is 
not the main concern in this work) is due to the non- 
linear interaction of Kelvin waves (helical perturbations 
of vortex filaments). It transfers kinetic energy down- 
scale, creating shorter and shorter waves, until, at very 
short length scales, phonon emission turns the kinetic en- 
ergy of the waves into the sound energy. (In ■^He-B, the 
dissipation in the zero-temperature limit is likely to be 
associated with a different, Caroli-Matricon mechanism 
of energy loss from Kelvin waves into the quasiparticle 
bound states.) In our calculation, a small numerical dis- 
sipation arising from the discretization and the vortex re- 
connection procedure plays the role of phonon emission 
[27t . If we continued the calculation for a longer time, 
this numerical dissipation would make the vortex length 
to decay, but the energy and the length can be considered 
effectively constant over the timescale of interest here. 

To highlight the presence of any structure in the vortex 
tangle shown in Figure [1] we convolve our discrete vortex 
filaments with a Gaussian kernel, and define a smoothed 
vorticity field u) as 

N s' 

'^(''' = E 7^37^ exp(-|s,: - r| V2a2)AC , (3) 
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where = dsi/d^ is the unit vector along a vortex at 
the discretization point = Si(^,t), and the smoothing 
length a is of the order of £. We test that, under this 
smoothing operation, a collection of randomly oriented 
vortex lines whose separation is of the order of i yields 
a; w 0; conversely, an organized bundle of vortex fila- 
ments of the same circulation yields a smooth vorticity 
distribution. Because of the wigglyness of the vortex lines 
at short lengthscales created by the Kelvin cascade, we 
also test the smoothing procedure on vortex lines after 
adding ten helical waves with an imposed k~'^^^ Kelvin- 
wave amplitude spectrum, as shown in Figure [3] 




^ from antiparallel (top) and parallel (bottom) vor- 
In the former, the average vorticity is small but 



FIG. 3: Volume rendering (a semi-transparent represen- 
tation) of the magnitude of the smoothed vorticity field 
|a;(r)|,s-^ 
tex pairs 

not zero due to the presence of small-amplitude Kelvin waves 
(note the wigglyness of the vortex filaments). In the latter, 
the contributions of the two vortex strands add up. 



Figure |4] shows the smoothed vorticity corresponding 
to Figure [H vortical 'worms', such as those seen in nu- 
merical simulations of ordinary turbulence, are clearly 
visible. 
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tions of quantum turbulence at relatively high temper- 
atures, 0.7K < T < Tx « 2.17 K, where thermal ex- 
citations form a viscous normal fluid component which 
interacts with the vortex lines via the mutual friction 
force. However, in this finite-temperature regime, vorti- 
cal structures belong to the normal turbulent fluid, and 
the friction force 36] would naturally induce superfluid 
bundles around them (although their stability is an open 
question). Here we focus on the development of vorti- 
cal bundles in the perfect superfluid, as a consequence 
of pure incompressible Euler dynamics. An intermit- 
tent vortex structure has also been noticed in a recent 
[33 | Gross-Pitaevskii simulation of turbulence in a Bose- 
Einstein condensate. However, in their calculation the 
vortex cores were (typically of a Bose-Einstein conden- 
sate, either realized numerically or in the laboratory) 
very close to each other: i/a^ sa 4.5 to 9, in contrast to 
our much larger £/ao ~ 3.5 x 10'^ typical of '^He and '^He- 
B superfluids; compressible effects (rapid density changes 
near vortex cores, sound waves, and the application of ar- 
tificial damping which affects vortex positions) which are 
absent in our (incompressible) Biot-Savart model were 
also likely to have played a role. 
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FIG. 5: The cross-section 2: = of the z-component of the 
smoothed vorticity field from Figure[4l ujz{r),s~^ . Open cir- 
cles show where vortex lines cross the plane from z < 0, and 
asterisks show those coming from 2: > 0. The clustering of 
vortex points of the same sign is evident, visible as regions 
of large positive and negative smoothed values of (shown 
blue and yellow). 



FIG. 4: Volume rendering of the magnitude of the smoothed 
vorticity field of Figure[Tl |u'(r)|, s~^. Note the vortical worm- 
like structures, similar to the coherent structures observed in 
ordinary turbulence. 



This is clear evidence of three-dimensional coherent 
superfluid vortex structures at T = which arise sim- 
ply by the fundamental equation of motion ([l} of vor- 
tex dynamics. There have been reports on the exis- 
tence of vortex bundles [13, HH in Biot-Savart simula- 



To verify that the structures in Figure |3] are indeed 
bundles of quantized vortices, we show in Figure [5] a 
two-dimensional cross-section of the z-component of the 
smoothed vorticity a; in the z — plane of Figure |4l 
We overlay the intersections of vortex filaments with the 
plane z = distinguishing the sign of ujz- It is appar- 
ent that the smoothed vortical structures shown in Fig- 
ure |4] are indeed small bundles of aligned vortex filaments, 
which appear as small clusters of positive and negative 
vorticity. There are typically 2 to 5 vortex points in each 
cluster. 
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It is important to check that the vortex bundles are 
physicahy distinct coherent structures of a weh-defined 
scale, rather than just a part of a purely random distri- 
bution, which would also contain structures of any scale 
which would become prominent after averaging. Cluster- 
ing in any two-dimensional system of points (or of any 
other discrete objects) can be confirmed and quantified 
using Riple y's JC- function [st} . or, alternatively, Besag's 
L-function [3g|. These functions, frequently used in ap- 
plied statistics, were also used to detect spatial correla- 
tions in two-dimensional physical systems (in particular 
in Monte-Carlo simulations and also in biological applica- 
tions of the Ising model). Ripley's /C- function is defined 
as 

^(^)= M ^ ^ m,<d), (4) 

where M is the number of points within the area A = 
(with D the size of the domain), dij is the distance be- 
tween points i and j , and the I{x) is unity if the condition 
X is satisfied and zero otherwise. However, the related 
Besag's function is more convenient for our purposes, de- 
fined as 

L{d) = ^K{d)/iT~d- (5) 

L = means complete spatial randomness, i < im- 
plies dispersion, and i > aggregation (clustering). We 
show in Figure |6] Besag's function for a large collection 
of two-dimensional cross-sections of Figure ID The re- 
sult confirms that vortex points of the same sign tend to 
cluster, hence the vortex bundles of Figure |4] represent 
distinct physical entities rather than merely an element 
of a purely random distribution of vorticity. 

We now proceed with the interpretation of our results, 
paying particular attention to the relation to the ordi- 
nary turbulence. According to our calculations, parallel 
vortex lines tend to come together locally, forming clus- 
ters of vortex filaments of the same circulation in any 
cross-section through the flow. If the vorticity distri- 
bution was two-dimensional, the well-known inverse en- 
strophy cascade of the two-dimensional turbulence would 
transfer vorticity to large scales, creating large-scale re- 
gions populated by vortices of the same sign. The two- 
dimensional inverse cascade is due to the conservation of 
the enstrophy, Q = Jy \uj\'^dV, the second invariant of 
the Euler equation in two dimensions. In three dimen- 
sions, only the first integral exists, the kinetic energy 
E = Jy i|vpc?y, whose conservation results in the di- 
rect cascade of the energy to smaller scales in both two- 
and three-dimensional turbulent fiows. 

To find whether there is a link between formation of 
three-dimensional coherent structures and tendency of 
vortices to cluster in two dimensions we performed addi- 
tional numerical experiments. The two-dimensional sys- 
tem of interacting point vortices in the inviscid fluid was 
introduced by Onsager [4^ in 1949, and is known as the 
Onsager point-vortex gas. This is where the concept of 
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FIG. 6: Besag's function, L{d), averaged over a random selec- 
tion of 1,000 two-dimensional cross-sections (parallel to one 
of the coordinate planes) of the vortex tangle of Figure |4] 
in the final, statistically steady state (black) and the initial 
state (red). Vertical bars show the one standard deviation 
range. Note that L increases with d sharply from L = at 
d = to its maximum value for d ~ 4 x 10~® m, and then 
varies smoothly with d, reaching the secondary maximum at 
d ~ 2 X 10"'^ m. The fact that in the final statel/(d) is sig- 
nificantly larger than zero, in the whole range of d explored, 
confirms the existence of coherent structures (clusters of vor- 
tex filaments in each cross-section) of various scales d, with 
some predominance of d « 2 x 10"'^ m. On the contrary, the 
initial state has L = within errors at all scales, consistently 
with the purely random nature of the initial conditions. 



inverse cascade was born. Onsager found that, statis- 
tically, the states of the two-dimensional system corre- 
sponding to the formation of large-scale coherent struc- 
tures (clusters of point vortices of the same polarity) for- 
mally have negative temperature (a detailed analysis was 
later provided by Montgomery [4l|). The ordering of 
vortices and the properties of the inverse cascade in the 
point-vortex gas have been extensively studied, with the 
numerical study of the formation of large-scale circula- 
tion patterns in the two-dimensional periodic domain by 
Montgomery and Joyce [i^l being one of the first studies 
of this kind. Recently some of the authors of the present 
report demonstrated [43| the clustering of point vortices 
of the same polarity in a small system of just a few hun- 
dred vortices, by computing the dis balance between the 
vortices and antivortices within subdomains of the pe- 
riodic two-dimensional computational domain. Here we 
have simulated the evolution of a much larger number 
(10,000) of point vortices in two dimensions. The sys- 
tem evolves from a statistically uniform state with zero 
net circulation where the number of vortices equals that 
of antivortices. Using Besag's function, we have con- 
firmed quantitatively the non-random clustering of vor- 
tex points. 

To check if coherent structures discussed above could 
be a fossil of the two-dimensional reverse cascade that 
survives in three dimensions, we explored a modified two- 
dimensional Onsager's system of point vortices similar 
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to that shown in Figure [5) We modelled the vortex 
motion out of the plane by removing randomly chosen 
vortex-antivortex pairs from the plane at a certain rate, 
and then randomly re-introducing them at random po- 
sitions elsewhere in the plane. For example, if a moving 
vortex loop in three dimensions crosses the plane, the 
number of point vortices in the plane first increases by 
two, and then is reduced by two as the loop leaves the 
plane. As these disturbances violate the local enstrophy 
conservation, they can be used to model, in a heuristic 
manner, local three-dimensional effects. We found that 
Besag's function in such a disturbed two-dimensional sys- 
tem still shows tendency to clustering typical of the in- 
verse enstrophy cascade, but weaker than in the pure two- 
dimensional case. The (positive) value of Besag's func- 
tion depends on the rate of vortex removal. It is therefore 
likely that the two-dimensional tendency to cluster is lim- 
ited by the fact that the enstrophy is not constant, but 
this tendency is not entirely suppressed entirely. 

When interpreting our results, it is also useful to 
compare them with recent calculations of Numasato et 
al. |39j . who observed only the direct energy cascade in 
a two-dimensional Bose-Einstein condensate. These au- 
thors solved the Gross-Pitaevskii equation; this model 
allows vortex reconnections, hence does not conserve the 
number of vortices (i.e., the enstrophy). These authors 
also argued that the Gross-Pitaevskii model is compress- 
ible (the condensate's density vanishes at the vortex axes 
and large-scale sound waves are easily excited, as noticed 
by White et al. j44|] ) , and the compressibility alone could 
explain the lack of the inverse cascade. Can compress- 
ibility be equally important in our case? As we already 
noted above, in a typical atomic Bose-Einstein conden- 
sate, either realized numerically or in the laboratory, the 
distance between the vortices, £, is only a few times larger 
than the vortex core radius oq, while ^/ao = O^IO"^) in 
the superfiuid helium. Therefore, compressibility is less 
likely to prevent the formation of coherent structures in 
the superfiuid turbulence. 

In the ordinary turbulence, the origin of coherent vor- 



ticity structures is often attributed to the roll-u p o f 
vortex sheets by the Kelvin-Helmholtz instability [45|. 
The superfiuid turbulence considered here has no vortex 
sheets, so the coherent structures observed here seem to 
involve a different mechanism. 

Our future work will explore if the number of vortex 
lines in the bundles increases as the circulation is reduced 
(with kinetic energy fixed), as one expects that the limit 
/i — should yield the classical behaviour where the 
range of vorticity values in coherent structures is wider 
than that observed in our numerical experiments. To 
better understand the dynamics of coherent structures, 
we also plan to investigate the time evolution of Ripley's 
and Besag's functions for the three-dimensional vortex 
tangle and to explore the time correlations of the fiow. 

In summary, our numerical experiments with very in- 
tense quantum turbulence at absolute zero have revealed 
that the vortex tangle contains coherent vortical struc- 
tures, or bundles of vortex lines, which arise from the 
Biot-Savart dynamics alone, and appear to be similar to 
the vorticity 'worms' observed in the ordinary turbulence. 
This result sheds new light on the relation between the 
ordinary and quantum turbulent flows, suggesting that 
their connection can be deeper than usually assumed. 
Given the relative simplicity of the quantum turbulence, 
this may provide new insights into the nature of turbu- 
lence. 
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